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Abstract 

Femtosecond pump-probe measurements find pronounced dimensionality de- 
pendence of the optical nonlinearity in cuprates. Although the coherent two- 
photon absorption (TPA) and linear absorption bands nearly overlap in both 
quasi-one and two-dimensional (ID and 2D) cuprates, the TPA coefficient is 
one order of magnitude smaller in 2D than in ID. Furthermore, picosecond 
recovery of optical transparency is observed in ID cuprates, while the recovery 
in 2D involves relaxation channels with a time scales of tens of picoseconds. 
The experimental results are interpreted within the two-band extended Hub- 
bard model. 

Transport and magnetic behavior of strongly correlated electron (SCE) systems, which 
depend on the lowest excitations, have been studied extensively in recent years. High energy 
electronic excitations, in contrast, have received less attention. Linear absorption studies, 
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performed for the copper oxide based SCE systems |]], do not reveal features associated 
with optically dark two-photon states. Femtosecond spectroscopic studies in quasi-two di- 
mensional (2D) cuprates |2||| have revealed the strong role played by magnetic excitations 
in ultrafast non-radiative relaxation processes, which is related to the large Heisenberg ex- 
changes in the cuprates . More recently, interest has shifted to nonlinear optical studies 
which provide information about both odd parity one-photon states and even parity two- 
photon states |9HiT|. 



An additional motivation to study nonlinear optical properties of SCE materials origi- 
nates from the continued quest for new materials with ultrafast and strong optical nonlin- 



earity. The mechanism of optical nonlinearity in SCE systems |T2| is different from that in 
band insulators, and it is conceivable that new SCE materials can be found with nonlinear 
optical properties that meet the requirements for technological applications. Initial support 
for this idea has come from the demonstrations of strong two-photon absorption (TPA) 
along with picosecond recovery of optical transparency in the quasi-one dimensional (ID) 
cuprate Sr 2 Cu03 ||, and of giant electro-reflectance in several ID Mott insulators [ItJ. An 
important new question that has arisen involves the role of dimensionality in SCE nonlinear 
optical materials. In conventional band insulators confinement of the Wannier exciton in ID 
can strongly enhance the optical nonlinearity [fl3 1 . Whether or not similar dimensionality 



effects occur in SCE systems has not been investigated systematically so far. In the present 
work we compare the results of the sub-picosecond pump-probe transmission measurements 
in ID and 2D cuprates. We observe significant dimensionality dependences in the magni- 
tude of the nonlinear response as well as the relaxation dynamics, which are explained using 
cluster calculations within the appropriate extended Hubbard Hamiltonian. 

To clarify the dimensionality dependence systematically, we selected three materials, 
which are well-known strongly correlated ID and 2D charge-transfer (CT) insulators. Specif- 
ically, Sr 2 Cu03 (see Fig. 1(a)) and SrCu02 (Fig. 1(b)) have single and weakly coupled double 
Cu-0 chains, respectively Since in SrCu0 2 the interchain coupling due to 90° Cu-O- 

Cu bonds is much weaker than the coupling along the 180° intrachain Cu-O-Cu bonds, this 



material can be also classified as a ID spin system |6|||. Sr 2 Cu02Cl2 (Fig. 1(c)) is a 2D 
cuprate with the same Cu0 2 network as in La2Cu04 with the apical oxygens of the latter 
replaced with CI ions Importantly, the lattice constants in the Cu-0 chain direction 
in these materials are almost same. All of these materials show optical gaps around 2 eV, 
which corresponds to CT excitation, and are nearly transparent below 1.5 eV (see Fig. 3). 
The strong antiferromagnetic exchange interaction between holes on neighboring Cu sites, 
J ~ 2000 - 3000 K for Sr 2 Cu0 3 and SrCu0 2 H and J ~ 1400 K for Sr 2 Cu0 2 Cl2 [§!, 
gives rise to wide spinon and magnon bands in ID and 2D materials, respectively. In spite of 
the one-dimensionality of Sr 2 Cu03, a spin-Peierls transition does not occur in this system. 

Single crystals of Sr 2 Cu03 and SrCu0 2 are grown by the traveling- solvent floating-zone 
method ||, while Sr 2 Cu0 2 Cl 2 is grown by cooling the stoichiometric melt |4]]. Thin flakes 
with thickness of 50-100 /im are cleaved out for transmission measurements. Two optical 
parametric generators pumped by a KHz regenerative amplifier generate pump and probe 
pulses. By using second harmonic and difference frequency generation techniques, the laser 
system can generate pulses with a temporal width of 0.2 ps and photon energy from 0.3 
to 1.5 eV. We measure differential transmission AT/T, where T is the transmission in the 
absence of the pump pulse, as a function of the pump-probe delay. The pump and probe 
beams are polarized along the Cu-0 chains (b axis in Sr 2 Cu03, c axis in SrCu0 2 and a axis 
in Sr 2 Cu0 2 Cl 2 ). 

The temporal behavior of the photoinduced absorption change at 290 K, AaL = - 
ln(l+AT/T), where L is sample thickness, are shown in Fig. 2 for all the three materials 
with the pump intensity I pump ~ 0.2 GW/cm 2 . The 2D material shows absorption change 
for both co- and cross-polarizations of the pump and probe beams, regardless of the angle 
between the polarization of the pump light and the crystallographic axes. We have found 
that AaL oc I pump up to I pump ~ 10 GW/cm 2 and I pump ~ 1 GW/cm 2 for ID and 2D 
materials, respectively. This indicates that the nonlinear effect is third-order in the light 
field up to these I pump . One can observe from Fig. 2(a) that the temporal profile of AaL in 
Sr 2 Cu03 consists of a prompt component, which is determined by the laser pulse duration, 
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and a slowly decaying component with characteristic time ~ 1 ps. From comparisons of the 
temporal profiles at various pump photon energies we conclude that the relative magnitude 
of the decay component decreases with decreasing pump photon energy. The temporal 
profile of AaL in SrCu0 2 (Fig. 2(b)) is similar to that in Sr 2 Cu03. On the other hand, in 
Sr 2 Cu0 2 Cl 2 , there exists a slower component with characteristic time ~ 30 ps as well as the 
faster component with ~ 1 ps (see the upper panel of Fig. 2). 

The prompt component of AaL is due to the coherent optical nonlinearity associated with 
TPA. In Fig. 3 we have shown the TPA spectra, plotted against 0J pump + uj pro b e , for all three 
materials. In all cases, one- and two-photon absorption maxima nearly coincide indicating 
that overlapping TPA and linear absorption bands is a common feature of cuprates. The 
TPA spectrum of Sr 2 Cu03 at 290 K is similar to that obtained at 10 K, with maximum 
f3 ~ 150 cm/GW at 2.1 eV 0. We further observe from Fig. 3 that while in the ID systems 
Sr 2 Cu03 and SrCu0 2 (3 is of the order of 100 cm/GW and the width of the TPA band is 
about 0.5 eV, in Sr 2 Cu0 2 Cl 2 the TPA coefficient is one order of magnitude smaller and the 
TPA band is broader. In particular, the low energy tail of the TPA in Sr 2 Cu0 2 Cl 2 continues 
down to uj pump + ijJ pro be ~ 1.5 eV, while the TPA coefficient in ID materials falls below 0.1 
cm/GW at 1.5 eV. 

In order to understand the differences in the TPA spectra between the ID and 2D 
systems, we adopted the two-band extended Hubbard model, which enables us to take into 
account the CT nature of the excited states explicitly, 

<ij,a i 

+ V J2 n i n 3 +J2 e i H i' i 1 ) 

<ij> i 

where c| CT creates a hole with spin a on site i, n ia = cl a Ci a , rii = J2a n ia, and < ij > implies 
nearest- neighbor sites. In Eq. ([]]), t is the hopping between Cu and O sites, [/, is the on- 
site Coulomb repulsion between two holes (different on Cu and O sites), V is the Coulomb 
repulsion between holes on neighboring Cu and O, and eo — ec« is the site energy difference 
between O and Cu sites. 
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It is instructive to first consider the qualitative difference between ID and 2D within 
Eq. (|T]). There exists a single mirror plane in ID, and optical transitions are between 
symmetry subspaces that are "plus" and "minus" with respect to this mirror plane. Consider 
now a three-atom segment OCuO. The single hole in the segment can occupy the Cu-site or 
either of the two O-sites. We denote these by the "cartoon" configurations 010, and 100 and 
001, respectively. The ground state, as well as the excited + symmetry two-photon state are 
superpositions of 010 and (100 + 001), while the - symmetry optical state is simply (100 
- 001). The dipole operator \x within Eq. ([!]) is X^r^ (we take electronic charge e = 1), 
where fl gives the vector location of each atom, and has nonzero matrix element between the 
configurations (100 + 001) and (100 — 001). The strength of the dipole coupling between 
the one-photon state and the ground (two-photon) state, /ioi (^12), then depends on the 
"overlap" between the (100 - 001) one-photon state and the (100 + 001) component of the 
ground (two-photon) state. 

Consider now the 2D case, where each plaquette consists of one Cu and four O-atoms, 
again with a single hole. There occur two mirror planes now, and the symmetry subspaces 
are + +, H — , — h and — with respect to these. Eigenstates in the + + subspace are now 
superpositions of the configuration with the hole on the Cu, and four configurations with 
the hole occupying the different O-atoms. This subspace contains both the ground state and 
a two-photon state. The optical states are in the H — and — h subspaces, are degenerate, and 
are still superpositions of only two configurations each: the H — eigenstate is a superposition 
of the two configurations with holes on the O-atoms to the left and to the right of the 
central Cu, while the — h eigenstate is composed of the configurations with the holes above 
and below the Cu. The dipole coupling between the one-photon and ++ states then can 
involve only two of the four configurations with holes on O-sites. Therefore, the transition 
dipole moments between the ground and one-photon states, /%., and one- and two-photon 
states, /Z12, are smaller in the 2D than in ID. Since a and (3 are proportional to yUgi an d 
/ i oi/ i i2! respectively |T3 |, we expect both linear absorption and TPA to be weaker in 2D 



than in ID, with the reduction in TPA strength larger. The above physical arguments for 



small-size cluster picture are valid even with a single electron, indicating that the weaker 
TPA in 2D can be explained without introducing charge-spin decoupling in ID, as has been 
previously claimed |T3| . 



We have verified the above qualitative reasonings by exact diagonalization studies of ID 
and 2D clusters containing 4 Cu atoms (see Fig. |4]). The 2D cluster chosen is the largest 
system that can be diagonalized exactly. The choice of the ID cluster size was based on 
the requirement that the number of holes are the same in ID and 2D. The parameters 
considered were \t\ = 1.4 eV, Uqu = 10 eV, Uq = 3 eV, V = 1 eV and eo — ec u = 2 eV. 
In Fig. |, we have shown the lowest energy levels, along with the mirror planes a in ID 
and <j x and a y in 2D. The very low energy excited states in the + subspace in ID and + + 
subspace in 2D in Fig. 4 are spin excitations that play no direct role in optical nonlinearity 
and are not discussed further. The calculated electronic structure and the identification of 
optically relevant states are consistent with the recent group-theoretical analysis based on 



the excitonic cluster model 11 



Our numerical results can be summarized as follows: (i) The dipole couplings of the 
ground state with the lowest pair of H — and — h states in 2D are one order of magnitude 
smaller than that with the next higher pair of H — and — h states in the 2D lattice. Accord- 
ingly, the optical states in 2D is a higher energy pair of H — and — h states (see Fig. 4). 
in 2D is 0.36 in our units, as compared to ~ 1 in the ID cluster, thereby explaining the 
weaker linear absorption in 2D (see Fig. 3). (ii) As seen in Fig. [|, in 2D, there exist even 
parity states, which are antisymmetric with respect to both a x and a y ( — subspace). TPA 
to these states should occur for cross-polarized pump and probe beams. We have observed 
the TPA with the cross-polarized configuration in our experiments and will discuss these 
elsewhere. Note that any interaction (electron-phonon interaction, intrinsic asymmetry due 
to crystal twinning) that leads to violation of this strict symmetry principle will cause weak 
TPA to — states even with co-polarized pump and probe. It is conceivable that the broad 
nature of the TPA in 2D in the high energy region and the persistence of weak TPA in 
the low energy region (~ 1.5 eV)in Sr 2 Cu0 2 Cl2 (see Fig. 3) are due to the - - states. As 
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seen in Fig. |], there also exists a — state that occurs below the true optical states. This 
state was not found in the previous one-band model |D| (presumably because the smaller 
lattice investigated there did not have the symmetry of the 2D lattice), (iii) In 2D, the 
most dominant two-photon states are the lowest and the third charge excitation states in 
the + + subspace, whose dipole couplings are 1.25 and 1.35 to the two one-photon states, 
compared to dipole couplings of 3.4 and 1.2 between the dominant two-photon states and 
the optical state in ID. The lower of the two + + excited states occurs slightly below the 
one-photon states in 2D for the parameters we have chosen, but with small modifications 
of these parameters it can also occur slightly above the one-photon states. Thus in both 
ID and 2D the dominant two-photon states are close to the one-photon states in energy, in 
agreement with experiment. The numerical simulation predicts TPA in ID larger by roughly 
one order of magnitude than in 2D even from these small cluster calculations. As discussed 
above, this is a consequence of the larger coordination number in 2D. (iv) We have included 
direct 0-0 hopping up to |t|/3 into Hamiltonian ([[]) and found no changes in our above 
conclusions. 

The other distinct dimensionality dependence is the difference in the relaxation time. 
Specifically, in Sr 2 Cu0 2 Cl 2 , there exist decay components with characteristic time ~ 30 ps 
and ~ 1 ps, while ID materials show only the latter component. The intensity of the decay 
components reduces with an decrease in pump photon energy (see Fig. 2). This indicates 
that the decay mechanism is associated with excitation of real carriers. We ascribe the 1 
ps relaxation to non-radiative channels through spinon (in ID) or magnon (in 2D) states, 
which exist below the optical gap and whose width is ~ 1 eV @||. We ascribe the slower 
30 ps optical relaxation in Sr2Cu02Ci2 to the existence of low energy H — and — h states 
with weak dipole coupling to the ground state and the forbidden - - electronic state that 
can act as trap states upon photoexcitation (see Fig. 4). Such subgap electronic excitations 
are absent in ID, and therefore the trapping of the optical excitation cannot occur. We 
note that the undoped 2D cuprates YF^CusOg and Nd2Cu04 have also shown two decay 
components with comparable decay rates @. 
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In conclusion, pronounced dimensionality dependence of the optical nonlinearity in the 
cuprates is found experimentally and theoretically. Coherent optical nonlinearity is domi- 
nated by two-photon states whose energy locations relative to the optical states are same 
in both ID and 2D materials. However, large nonlinearity and ultrafast relaxation are 
characteristics of only ID SCE systems, which make them promising materials for ultrafast 
optoelectronics |§. The spin excitations with large energies, promoting the relaxation, are 
unique to SCE systems. The smaller nonlinearity in 2D is due to the larger coordination 
number of the atoms, and the occurrence of a slower relaxation channel in 2D is most proba- 
bly associated with a fundamental difference in the electronic structures in ID and 2D, viz., 
the occurrence of electronic energy states below the optical states in 2D, and their absence 
in ID. 
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FIGURES 

FIG. 1. Crystal structures of cuprates. (a) Sr 2 Cu0 3 (b=3.91 A), (b) SrCu0 2 (c=3.92 A), (c) 
Sr2Cu02Cl2 (a= b=3.98 A). Schematics of the Cu-0 networks are also shown. 

FIG. 2. Temporal profiles of AaL, which are normalized at their maximum values, at 290K in 
Sr 2 Cu0 3 (a), SrCu0 2 (b), and Sr 2 Cu0 2 Cl 2 (c) for (pump + probe) energies of (1.46 + 0.95 eV) in 
the upper part and (0.95 + 1.03eV) in the lower part. One can observe that the picosecond decay 
components are pronounced for the higher pump photon energy. The upper panel shows AaL, now 
plotted on a logarithmic scale, for the upper cases ((a), blue; (b), green; (c), red) plotted against 
extended time scale. The black solid line shows a fitting curve to the AaL of Sr 2 Cu0 2 Cl 2 with 
two exponential functions with decay times of r=l ps and r=30 ps. 

FIG. 3. TPA coefficient (3 versus 0J pump + Lo pro be at 290K in Sr 2 Cu03 (a), SrCu0 2 (b) and 
Sr 2 Cu0 2 Cl 2 (c). Circles, squares and triangles correspond to pump energies at 1.46, 1.31 and 0.95 
eV, respectively. Solid lines indicate the linear absorption a. 

FIG. 4. Cluster models studied numerically and the corresponding energy levels in ID (a) and 
2D (b). Filled (open) circles are Cu (O) sites. Dashed lines denote mirror-plane symmetries used. 
States labeled * (1-photon states) have large dipole coupling to the ground state, those labeled # 
(2-photon states) have large dipole coupling to 1-photon states. See text for parameters. 
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Fig. 1. M. Ashida et al. 
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